Tkachenko modes as sources of quasiperiodic pulsar spin variations 
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We study the long wavelength shear modes (Tkachenko waves) of triangular lattices of singly 
quantized vortices in neutron star interiors taking into account the mutual friction between the 
superfluid and the normal fluid as well as the shear viscosity of the normal fluid. The set of 
Tkachenko modes that propagate in the plane orthogonal to the spin vector are weakly damped 
if the coupling between the superfluid and the normal fluid is small. In strong coupling, their 
oscillation frequencies are lower and are undamped for small and moderate shear viscosities. The 
periods of these modes are consistent with the observed ~ 100 — 1000 day variations of spin for PSR 
1828-11. 



I. INTRODUCTION 



During the last decade mounting evidence has emerged 
for the existence of long-period oscillations in a handful 
of pulsars [l|, H, HI, EL HI- An outstanding example of 
this kind of phenomenon is observed in PSR 1828-11 [|[. 
Its timing residuals are modulated with periods of 256 
and 511 days, while a 1009-day periodicity is inferred 
with lower confidence. These timing residuals coincide 
with periodic modulations of the pulse shape, which is 
a strong indication for the precessional motion(s) of the 
pulsar with the periods quoted above In fact, var- 
ious models of precessing neutron stars @, 0, H, m, E3| 
fit the timing data fairly well. However, it is known 
that free precession is incompatible with the existence 
of a superfluid in the pulsar's interior if the superfluid 
is strongly coupled to the normal fluid [ill G2, 0> E3 • 
The mutual friction for the superfluids present in the 
crust and core of neutron stars, which are derived us- 
ing microscopic calculations, cover a broad range of val- 
ues. A reliable calculation necessarily involves the super- 
conducting and superfluid properties of the fluid(s) and 
the non-superconducting material, which resists to vor- 
tex motion. While precession remains a viable model for 
the quasiperiodic oscillations observed in pulsar timing 
data, here we follow a different route [l5| by exploring 
the propagation of Tkachenko modes in pulsars as the 
source of long term variations. 

Charge neutral superfluids in neutron star interiors ro- 
tate by forming an array of singly quantized vortices. 
In their lowest energy state the vortices form a two- 
dimensional triangular lattice. The lattice supports col- 
lective elastic modes, Tkachenko waves, in which the vor- 
tices are displaced parallel to each other [l6|, [l?], [H, [l|| . 
Their undamped propagation would lead to variations of 
the angular momentum of the superfluid due to the lo- 
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cal variations of the density of the vortex lines as well as 
periodic variations in the rotation and spin-down speeds 
of the star. 

Ruderman [2(| pointed out that the frequency of the 
Tkachenko modes is of the order of several hundred days 
and these modes could be responsible for the quasiperi- 
odic timing residuals observed in the Crab pulsar. Lit- 
tle attention has been paid to the role played by the 
Tkachenko modes in neutron stars since Ruderman's 
1970 paper. Here we reinvestigate the propagation of 
Tkachenko modes in neutron star superfluids and calcu- 
late how these modes are damped by mutual friction and 
the shear viscosity of normal matter. These factors are 
clearly important for the continuous propagation of these 
modes, which would certainly lead to observable effects. 
Provided that the observed variations are indeed caused 
by the shear modes of the lattice we can ask the ques- 
tion: what do the quasi-sinusoidal variations tell us about 
the microscopic physics governing superfluids in neutron 
stars? 

This question can be answered by describing the 
physics present at micro- and mesoscopic scales (the up- 
per limit on these scales is set by the size of the neutron 
vortex) in terms of the few parameters that enter the 
equations of superfluid hydrodynamics. These parame- 
ters include the kinetic coefficients and other local char- 
acteristics (e. g., density, baryon and lepton fractions, 
etc). The predictions obtained within this hydrodynamic 
description can be presented in a form that is indepen- 
dent of the details of physics at the micro- and mesoscopic 
scales. 

In this paper we pursue this "model independent" 
approach, although some comments about the physics 
present at intermediate scales are in order. A key prob- 
lem concerning the crust of neutron stars is the question 
of whether neutron vortices are pinned to the crustal nu- 
clei both under static and dynamical conditions. The 
static problem always has pinning solutions for infinites- 
imally attractive interactions between a neutron vortex 
and a nucleus (for repulsive interactions the vortex is 
localized in between the lattice sites) [2l|, [22| . The dy- 
namical problem of repinning under the action of exter- 
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nal superflow [23[ allows for repinning solutions only for 
strong pinning potentials and current estimates for the 
pinning energy do not favor repinning solutions. How- 
ever, if the axisymmetry of the problem is lost, as in the 
case of processing neutron stars, the vortices may not pin 
at all [H. 

Another issue concerning neutron star crusts is that 
they may often undergo starquakes [25l . |26( , which drive 
the vortex lattice out of its equilibrium state. In fact, the 
form of the lattice may change from a triangular shape 
into some other structure reflecting the symmetries of 
the underlying nuclear lattice. Moreover, other possible 
scenarios include cases where the vortices may acquire 
multiple quanta of circulation or the lattice may contain 
impurities and vacancies. Thus, our working assumption 
that the ground state of the lattice consists of a simple 
triangular array of vortices may not be applicable to the 
crusts of frequently quaking neutron stars. 

Regarding the cores of neutron stars, a majo r ques- 
tion is whether the proton fluid is a type-I I27I. 28 
or type-II superconductor [3(| HH [H, [H, l34l. l35l. 136 
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37. 



l38l . l39i . |40( (it may also become unpaired at high densi- 
ties) . The friction in the case of type-II superconductors 
is large (38|, which implies that precession is impossible. 
If, however, the proton fluid is a type-I superconductor 
the friction is found to be small and compatible with 
precession |15| . 

Finally, we would like to point our that the superflu- 
ids may develop quantum turbulence when the velocity 
difference between the superfluid and the normal fluid 
exceeds a certain critical value. Turbulent states in su- 
perfluid He may be generated in many ways, e. g., by 
driving the normal fluid along the vortex lines through 
a temperature gradient. Neutron stars superfluids may 
also develop turbulent states [4l|, l43l. which could 
be caused either by unstable precession [4J| or tectonic 
activity in the crusts |45| . 

This paper is organized as follows. In Sec. |TT] we 
recapitulate the equations of superfluid hydrodynamics 
that include the combined effects of vortex tension, mu- 
tual friction, and shear viscosity. Sec. IIIII is devoted 
to the derivation of the characteristic equation for the 
Tkachenko and inertial modes and their numerical study. 
Our conclusions are summarized in Sec. HVl 



II. 



SUPERFLUID HYDRODYNAMICS WITH 
VORTEX TENSION 



It is convenient at this point to discuss the three main 
length scales that appear in our study in more detail. 
The relevant length scales are the vortex core radius 
~ 10~ 12 cm, the intervortex spacing ~ 10~ 3 cm, and 
the size of the superfluid phase ~ 10 5 cm (we refer to 
these three different scales as micro, meso, and macro- 
scopic scales). A hydrodynamic description requires av- 
eraging over the mesoscopic scales. In order to describe 
the deformations of the vortex lattice, one needs an ad- 



ditional dynamical variable, the local deformation of the 
vortex lattice e(r), which can be incorporated into the 
Bekarevich-Khalatnikov superfluid hydrodynamics [4t| . 
A hydrodynamic description of superfluids that includes 
lattice deformations has been studied by a number of 
authors [13, EH, H^, [5(| and in this work we use the 
Baym-Chandler version of superfluid hydrodynamics [49l | 
to study the Tkachenko modes and their damping due to 
mutual friction and shear viscosity. 

The fluid motions are naturally separated into the cen- 
ter of mass motion (known as first sound) and second 
sound, which corresponds to temperature dependent rel- 
ative oscillations between the superfluid and the normal 
component. These motions are conveniently described in 
terms of the total mass current j = pn^n + Ps^s and 
the relative velocity w = v^v — vs, where ps (pn) and 
vs (vat) are the superfluid (normal fluid) density and ve- 
locity, respectively. In addition, one needs an equation 
for the time variations of the lattice deformation e(r). 
We note in passing that p = p^ + ps is the total mass 
density and the tiny effects arising due to the inertia of 
the vortex lines are neglected. 

The linearized version of the fundamental superfluid 
hydrodynamic equations written for the net mass cur- 
rent, the relative velocity, and the superfluid velocity are 
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where ft = (0, 0, f2) is the spin vector (f2 is the pulsar 
rotation frequency), P = Pq — p(£l x r) 2 /2, Pq is the 
pressure in the fluid at rest, and a is the vortex elastic 
force density defined as 



a- = ps [2V± ■ (V± ■ e) 



V 2 ± e] 



2QX 



dz 2 ' 



(4) 



where Vj_ is the gradient in the x — y plane and ps — 
p5fi.fi/8m7v is the shear modulus of the trian gula r vortex 
lattice calculated by Tkachenko 0, [13, El 1H ■ More- 
over, mj\r is the bare neutron mass and the vortex tension 
is given by 



(5) 



where a is the coherence length and b = {irhj 'v^mArfi) 
is the vortex radius of the triangular lattice. The New- 
tonian gravitational potential <fi satisfies the equation 



V 2 4> = V z 



+ 4> N ) = 4:nG(ps + Pn), 



(6) 



where G is the Newton's constant and <ps and <pN are the 
gravitational potentials of the superfluid and the normal 
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fluid, respectively. The force density C is defined as Ci 
VfcTjfc, where is the viscous stress tensor 



-77 ( ViV Nk + VkVNi - -^S ik V • v w ) , (7) 



whereas r\ is the shear viscosity. Note that we do not take 
into account the effects from bulk viscosity and thermal 
conductivity. Finally, the mutual friction force is 



n x 



u> x 



de 



+P'PS 



u> x 



de 



(8) 



where u> — V x vg is the quantum circulation vector, 
n = ijj/u), and and (3' are the phenomenological mutual 
friction coefficients. 

The dissipative terms in the hydrodynamic equations 
defined above (such as the mutual friction forces) are 
the most general expressions that can be used in the 
description of a rectilinear vortex lattice in equilibrium 
that are still compatible with the conservation laws and 
the assumption that the dissipative function is a posi- 
tive quadratic form of the perturbations (higher order 
terms are neglected). Furthermore, we would like to re- 
mark that Baym-Chandler hydrodynamics only describes 
the linear order corrections in the lattice displacements, 
which are considered to be small. The neutron vortex 
lattice in neutron stars may, however, become unstable 
towards forming a tangle with turbulent superfluid flow, 
which would require a thorough revision of the hydrody- 
namic equations, fn particular, the form of the mutual 
friction force © should be changed to the one suggested 
by Groter and Mellink [51| where f cx w. These aspects 



of superfluid dynamics are beyond the scope of this work 
(see Refs. [M \M, E|). 



III. OSCILLATION MODES 

We consider plane wave perturbations with respect to 
the equilibrium, which corresponds to uniform rotation. 
We use a Cartesian system of coordinates where the z- 
axis is directed along the spin vector u>. The vectors j and 
w can be decomposed into transverse and longitudinal 
parts, i.e., j = j t + j; and w = w t + w/. The transverse 
parts we are interested in satisfy the condition 



V ■ jt = V ■ w t = 0. 



(9) 



The perturbation equation for the transverse components 
of the vectors j and w derived from Eqs. CD)-© are 
(hereafter the subscript t is suppressed) 



)Pa = 0, (10) 



dji 
dt 

+ (2ei mn Q m w n - — - ft) Pa = 0, (11) 
at ps 

Ov' d€ 

-^+2e lmn n m -^-P a = 0. (12) 

In the equations above we used the projector Pa = 
Sn — kiki/k 2 where k is the wave vector. Our coordi- 
nate system is such that the wave vector lies on the z — x 
plane, i.e., k = (A: sin 9, 0, A; cos ff) where 6 is the angle 
formed by the vectors fl and k. 



Writing the time perturbations as ji(t) 



Ji 



(a similar definition is used for the other vectors) we 
obtain, after some algebra, the characteristic equation 



det \\K, 



where 
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isiNg -f)ls 
■P*h p + P(d + -y N g) -P*(l- lN h) 
Ph P*(d + lN g) p + P{l- 7N h)J 



(13) 



We used the following shorthand notations in the defini- 
tion of the matrix Kiji 7jv/s = Pn/s/P: d 1 ^ 2 — cos 9, 
(3* = 1 -P',fj = vk 2 /(2Qp), a = (4 - d)/3, (3 = 7^/?, 
and $' = 7^ 1 /3 / . Moreover, we defined 

9 = 77§— i^s ~ d(ps - 20A)] , (14) 
h = Sr~ ^ " d ^s + 2QX)] . (15) 



the density because ps ~ Ps an d X ^ ps- The den- 
sity appears only through the normalization of the shear 
viscosity f] ~ P ■ The eigenmodes of the matrix in 
Eq. (fl3|) provide the oscillatory modes in the general case 
where the shear viscosity of normal matter and the mu- 
tual friction are included. In the non-dissipative limit 
(/3 = j3' = T) = 0) the modes separate into two indepen- 
dent sets that describe the inertial and Tkachenko modes, 
respectively. The (real) eigenfrequencies of these modes 
in units of 2S7 are 



Note that the coefficients g and h are independent of 



PI = ±id^ 2 , PT = ±i{(d + g)(l-h)} 1/2 , (16) 
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FIG. 1: (Color online) Dependence of the period P — 2tt/\pt\ 
of the Tkachenko modes on the wave vector for d = (solid 
black line), d = 10 -12 (dashed red line), and d — 10~ 6 
(dashed-dotted blue line). For large wavelengths the periods 
are on the order of 100 days. In particular, for SI* = 15.51 
Hz we obtain that P(k m i n , d = 0) = 331 days if the core size 
is R c — 10 km and 256 days when R c = 7.7 km. 



where the indices / and T refer to inertial and Tkachenko 
modes, respectively. 

If the Tkachenko modes are generated within super- 
fluid shells with the width of R c ~ 10 km their corre- 
sponding wave vectors are of the order of k mm = 2ir/R c ~ 
6.28 x 10 -6 cm, which sets the lower limit on the wave 
vector. Since the hydrodynamic description breaks down 
at length scales ~ 106 (b ~ 10~ 3 cm is the intervortex 
distance) the wave vector is bounded from above by the 
value fcmax = 73.3 cm -1 . We are interested in the small 
wave vector limit k ~ /c m j n that describes vortex density 
waves across the entire superfluid shell. 

The parameters g and h are of the order of s = 
(h 2 k^ [lin /2m n )(8h£l)~ 1 . For instance, when k = /c m i n and 
ft* = 15.51 Hz (the rotation frequency of PSR 1828-11) 
we obtain that s — 10~ 14 . Therefore, for s <C 1 the eigen- 
modes corresponding to Tkachenko waves in the dissipa- 

1/2 

tionless limit are given by p T = ±i [(d + g)] 1 . In the 
limit d <C g where the wave vectors are highly collinear 
to the spin vector we obtain that px = Hy/g and, in 
the opposite limit d 3> g, the Tkachenko modes become 
identical to the inertial modes, i.e., Pt = Pi- 

Fig. Q] displays the period of the Tkachenko modes 
without dissipation as a function of their wave vector. 
Only the long wavelength perturbations have periods on 
the order of 100 days, which are then relevant for obser- 
vations. In this limit the periods rapidly decrease for per- 
turbations with finite d. The period P(k min , d = 0) = 331 
days for tt+ = 15.51 Hz suggests that the shortest of the 
periods observed in PSR 1828-11, which corresponds to 



256 days, should be identified with the fundamental os- 
cillation mode. Oscillations with larger periods should 
then be identified with the higher-order harmonics of this 
mode. A period of 256 days can be obtained by adopt- 
ing R c = 7.7 km, which translates into k = 8.16 x 10~ 6 
cm -1 . This value is close to the upper limit for the size 
of any superfluid region inside a medium-heavy neutron 
star. It is important to remark that if the width of the 
superfluid region where Tkachenko waves can be found 
is at least one order of magnitude (or more) smaller than 
the values of R c quoted above the modes will be too fast 
to account for the long-period oscillations observed in 
PSR 1828-11 (as long as the renormalization effects due 
to mutual friction discussed below are neglected). 

We now consider the effects of the shear viscosity and 
the mutual friction on the propagation of Tkachenko 
modes. It is convenient to use the drag-to-lift ratios £ 
and (J instead of (3 and (3' to describe mutual friction. 
These ratios are related by the following equations 

= i ^i- ^-co (17) 

Microscopic calculations indicate that £' ~ 0. The limit 
C — * corresponds to weak coupling between the vor- 
tices and the normal fluid, while C — > o° implies strong 
coupling. 

Figure [5] shows the dependence of the modes derived 
from Eq. (| 13|) on the drag-to-lift ratio £ for several values 
of the shear viscosity and d = 0. The value of f\ is deter- 
mined assuming a constant density of 3 x 10 14 g cm' 3 . In 
the limit where C and rj vanish we recover the results for 
the non-dissipative case discussed above. For 77 = the 
real part of the Tkachenko mode, which is doubly degen- 
erate, vanishes only in a narrow window of values of £. 
For larger values of C, which corresponds to the strongly 
coupled region, the real part reaches an asymptotic value 
that is about 25% smaller than its value in the undamped 
limit. Note that in our plots only the regions where the 
modes change significantly are shown. 

Assuming that the normal fluid is inviscid, the results 
in Fig. pimply that there are oscillations with even longer 
periods in the strongly coupled limit. The Tkachenko 
modes are significantly damped by mutual friction in the 
region where the real part vanishes. There the number of 
imaginary roots of the characteristic equation increases 
by one. Moreover, one of the imaginary roots is given 
by Im u) = i/30,, which continues beyond the figure's y 
scale. This reflects the damping of the differential rota- 
tion between the superfluid and the normal fluid caused 
by mutual friction. This damping has no effect on the 
Tkachenko modes in strong coupling. For moderate val- 
ues of viscosity (77 = 5 x 10 17 dyn s cm~ 2 ) the real part of 
the Tkachenko mode is reduced in the strongly coupled 
region. In this case, its imaginary part is smaller than 
the real part and, therefore, the oscillations are weakly 
damped. For large values of the shear viscosity (rj ~ 10 19 
dyn s cm~ 2 ) the real part of the Tkachenko mode van- 
ishes in the strongly coupled limit. Finally, note that 
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FIG. 2: (Color online) Dependence of the real (solid black 
line) and imaginary (dashed red line) parts of the Tkachenko 
modes (uj = ip) on the drag-to-lift ratio Q for r\ = (upper 
panel), r\ = 7.5 x 10 17 (middle panel), and r\ = 1 x 10 19 (lower 
panel) in dyn s cm" 2 units. All modes are normalized by the 
non-dissipative value of the Tkachenko mode Qt = 2n/P. 
The modes are computed taking d = 0, which means that 
there are no inertial modes. 
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FIG. 3: (Color online) The same as in Fig. [2] in case of d = 
10" 14 and 77 = (upper panel), 77 = 10 19 (middle panel), and 
77 = 5 x 10 19 (lower panel) in dyn s cm -2 units. Note that 
both the Tkachenko and inertial modes are displayed here. 



larger because 77 ~ T 



IV. CONCLUSIONS 



there are no inertial modes when d = 0. 

The modes when d — 10" 14 are shown in Fig. [31 As dis- 
cussed above, in the non-dissipative limit the Tkachenko 
and inertial modes coincide for sufficiently large d. For 
77 = the modes can be distinguished in the strongly 
coupled limit because the Tkachenko mode vanishes for 
sufficiently large values of (. When larger viscosities are 
considered (77 > 10 19 dyn s cm" 2 ) the difference between 
the real parts of the inertial and Tkachenko modes can 
be clearly resolved. If we increase 77 even further we see 
that the real part of the inertial mode decreases and the 
imaginary part, which increases with 77, becomes rele- 
vant. Finally, the real part of the inertial mode vanishes 
at 77 ~ 5 x 10 19 dyn s cm" 2 . 

The outer cores of neutron stars are mainly composed 
of light baryons, which pair in the isospin triplet states, 
and leptons. For densities of 2 — 3 x 10 14 g cm -3 and 
temperatures of T ~ 10 8 K the shear viscosity of the 
electron fluid was determined to be in the interval be- 
tween 8 — 40 x 10 17 dyn s cm" 2 52]. This value of the 
temperature is a realistic upper bound on the tempera- 
ture in the core of neutron stars except for very young 
objects such as the Vela and Crab pulsars. For colder 
stars the viscosity could be a few orders of magnitude 



Our results indicate that Tkachenko modes are broadly 
consistent with the weakly coupled theories between the 
superfluid and the normal fluid, independent of the value 
the shear viscosity. The subclass with d = cos 2 9 = 
has periods that are consistent with the lowest observed 
periodicity in PSR 1828-11 of 256 days. 

The existence of Tkachenko modes in the strongly cou- 
pled region depends on the shear viscosity of normal mat- 
ter. For low viscosities the Tkachenko modes are (in 
strong coupling) renormalized to values that are a few 
times smaller than their non-dissipative limits. This im- 
plies that in strong coupling the Tkachenko oscillations 
have periods that are larger than their non-dissipative 
counterparts. In fact, the damping caused by mutual 
friction is not always strong enough to preclude an os- 
cillatory behavior. Therefore, we conclude that the long 
term variation in the spin of PSR 1828-11 can in princi- 
ple be explained in terms of Tkachenko oscillations within 
superfluid shells for a broad range of values of the mutual 
friction and the normal fluid shear viscosity. Finally, we 
note that larger wave vector oscillations corresponding 
to periodic motions on shorter length scales may lead to 
phenomena that are responsible for the observed timing 
noise in pulsars. 

Our model necessarily involves certain approximations. 
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For instance, we have adopted the two-fluid superfluid ACKNOWLEDGMENTS 
hydrodynamics, which should be modified in order to 
account for the m ultip le fluids in the neutron star's 
core [H, [H [H HI SaH, SI . Furthermore, the cylin- 
drical symmetry of our setup and the assumption of uni- 
form density need to be relaxed in more realistic treat- 
ments of spherical superfluid shells with density gradi- J.N. acknowledges support by the Frankfurt Interna- 
ents. tional Graduate School for Science (FIGSS). 
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